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Abstract 

I give a brief review of results obtained recently at Ecole Normale 
on the depinning transition of interfaces and contact lines, using a 
variety of approaches: non-local Monte Carlo algorithms, dynamical 
renormalisation group calculations to 2-loop order, and exact solution 
of an infinite-range model. 
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1 Introduction 

The motion of interfaces in random media (or more generahy of " elastic man- 
ifolds"), and in particular the existence of a depinning threshold analogous 
to a critical point, is a now classic problem of statistical mechanics which ap- 
pears under various forms in many areas of condensed matter physics Q . In 
spite of much progress during the last 15 years [|| many interesting questions 
remain open and I take this opportunity to present some results recently ob- 
tained by several groups at Ecole Normale. A personal motivation to talk 
about this topic here is that among many domains of common interest with 
Gene Stanley this is one of the few where I can claim priority, having studied 
contact lines on disordered substrates with Y. Pomeau Q before Gene's own 
work on interface depinning Q - see ||5| for relevant references up to 1995. 

*Laboratoire associe au CNRS ct aux Universites Paris VI et Paris VII. 
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The incentive to look at that problem again came initially from experi- 
mental results obtained in our laboratory by E. Rolley and his group on the 
helium-cesium contact line 0]- A main goal was to improve the deter- 
mination of the depinning threshold in order to obtain more precise values 
for the roughness and velocity exponents, in the hope to reconcile numerical 
simulations with the experimental results and with new theoretical predic- 
tions due to LeDoussal and coworkers As will become apparent from 
this brief review the results brought some surprises and call for further work. 



2 Algorithms and simulations 

A detailed description of the motion of a fluid-solid contact line is quite 
intricate if one tries to take into account in a realistic way both the hydro- 
dynamic aspects @ and the substrate inhomogeneity. For this problem and 
for interface depinning in general most theoretical studies therefore focus 
on the effects of disorder close to the threshold, when the mean velocity is 
sufficiently small for inertial effects to be negligible, at least as a first approx- 
imation. It is also customary to disregard the possibility for the interface to 
develop overhangs in order to pass round unfavorable regions. With these 
approximations the equation of motion at zero temperature may be written 



du/dt = K[u] + Fe^t + r/(n) , (1) 



where u{r, t) is the displacement measured from a flat reference (d — 1)- 
dimensional interface. f] is an elastic force that tends to keep the 

interface flat, its form depends on the system considered and is non-local for 
the contact line. Ff^^t is an external driving force and r/(n) represents the 
interaction with the disordered medium or substrate, it will be restricted in 
the following to random, spatially uncorrelated pinning forces distributed 
according to a given probability density p{r]). 

An essential difference with the much studied Edwards- Wilkinson and 
Kardar-Parisi-Zhang growth models Q is that here the disorder is quenched, 
as the random part of the force felt by an immobile interface element re- 
mains constant in time. This non-linearity makes the problem quite difficult 
to tackle analytically. Numerical simulations of (|l]) with continuous space 



and time also raise technical difficulties |1C]. In particular the interface will 



always get pinned by rare strong local defects if the maximum distorsion 
\ui — Uj\ allowed in a local move is fixed and if the random force is drawn 
from a continuous distribution, e.g., a Gaussian as used in field-theoretical 
calculations. Most simulations have in fact been performed on discretized 
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cellular-automata type versions of the problem, where the transverse posi- 
tions Ui{fi,t) of interface sites are updated at regular time intervals, the 
jumps being continuous or discrete. Such simulations are rather straightfor- 
ward to perform, at least in low dimensions, but even for these simplified 
versions it is not easy to obtain accurate values of the depinning threshold 



and of the associated critical exponents [11, 12 1. The interface velocity be- 
comes extremely small and fluctuates strongly, so very long simulation times 
and large system sizes are necessary. Another problem is that it is difficult 
to identify the equation that results in the continuum limit. 

This led Rosso and Krauth |13| to consider non-local Monte Carlo al- 



gorithms that bypass the dynamics and give direct access to the stationary 
configurations of the continuum equation (|l|), in the pinned phase. When 
the external force is increased the last configuration to get destabilized is 
the critical interface, it is unique for a finite sample with perodic boundary 
conditions. Rosso and Krauth showed how this interface may be efficiently 
determined using an algorithm in which one looks at each time step for the 
"front" of minimal length that it is favourable to move by one lattice unit. 
The roughness exponent C, for the interface width w as a function of its 
lateral extension, w ~ may thus be obtained right at the critical point 
with better accuracy than through the direct dynamical approach. 

For short-range elastic interactions with bounded distorsions the result 
is C — 0.63 in d = 2 in very good agreement with the value obtained 



from simulations on cellular-automata models [^, 15 1 which had led to the 
conjecture 

^ = zy_L/zy|| = 0.633 ±0.001, (2) 

where and are the transverse and longitudinal exponents for directed 
percolation. On the other hand, if the local distorsions are unconstrained 
and the restoring force is linear (i.e., the elastic energy is purely quadratic 
with respect to these distorsions), giving the so-called "quenched Edwards- 
Wilkinson" model, the exponent is found to be Qqew = 1-17, in agreement 
with other earlier studies |10|. 



If the restoring force is non-linear, as obtained for instance by adding a 
quartic term of the form (tij — Uj)^ in the elastic energy of the QEW model, 
one recovers C, ~ 0.63. This value appears to be robust and to correspond to 
a broad class of systems. The surprise is that it is obtained in the absence 



of an anisotropy term [17|, while the "standard" value C ~ 1.2 seems to 
be specific to the harmonic QEW model. It remains to take into account 
such higher-order terms in the RG approach and to determine precisely the 
different universality classes and their extent. 
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Figure 1: Critical interfaces obtained using the "variant Monte Carlo" al- 
gorithm, for short-range interactions with restricted height differences (top) 
and for long-range interactions (bottom). The size is = 512 in both cases 
(courtesy of A. Rosso). 



For the contact line problem the capillary tension acting on the fluid 
surface leads to an effective long-range interaction between the interface 
elements Q and the elastic force has the form 

. /■ , / u{r,t) -u{r',t) 
^ J dr ^^—^2 • (3) 

The simulations are more time-consuming in this case and a specific algo- 
rithm has to be devised, but the critical interface may still be determined 
accurately |18]. It is found to be much less rough, see figure (||), with 



Clr = 0.388 ± 0.002. (4) 

This result is also surprising, as experiments on well-controlled strongly 
disordered substrates Q give Cexp = 0.56 it 0.03, much higher than initially 
found on weakly disordered but ill-controlled substrates. The calculation 
of C using dynamical-RG methods is tricky, and the most recent work has 
shown that contributions at 2-loop order do not vanish |^] , contrary to initial 
expectations. One finds C^^^ = 0.47 at this order, leading to an extrapolated 
estimate (rg = 0.5 it 0.1, which lies closer to the experimental results but 
further from the numerical ones than the one- loop value C^^^ = 1/3 H. 
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3 An exactly solvable model 



In the case of equilibrium phase transitions the critical exponents of systems 
with power-law interactions such as (^) are intermediate between those of 
short-range systems and the mean-field results, which are obtained in the 
limit of weak uniform infinite-range interactions. The latter are also used as 
the starting point for field-theory RG treatments and it is natural to study 
analogous models for moving interfaces [Q, |2^. For discretized models the 
elastic restoring force on site i is simply given in the infinite-range limit by 
a sum over all sites 

Ki[u] Qc ^ {uj-Ui)= u-Ui, (5) 

j=l,N 

SO the motion of an interface element is just that of a particle in a one- 
dimensional potential, but with an additional coupling to the position u{t) 
of the center of mass, which is to be determined self-consistently. 

It is usually assumed in order to obtain analytical solutions that the 
infinite-range interactions have the effect of averaging out the disorder- 
induced fluctuations. The interface velocity v = du/dt is then independent 
of time in the mobile phase, in the limit of an inflnite system, and obeys a 
self-consistent equation. The solution shows that 

vr^ {Fe^t-F,)\ (6) 

where is the critical force needed to depin the interface. For non-singular 
distributions of the random pinning forces the critical velocity exponent is 
found to be H 

Omf = 1. (7) 

One would like to have a rigorous proof of this simple and physically rea- 
sonable result, but this is not so easy to obtain, as it involves an interchange 
of limits when time and size go to oo which may be tricky. Indeed, we have 
been able to solve exactly a particular model in the inflnite-range limit [^] 
and the exponent found in that case differs from the result 

3.1 The infinite-range Leschhorn model 

This solvable model belongs to a family of cellular-automaton type models 
due to Leschhorn in which the interface sites may only move forward 
by a discrete distance at regularly spaced time intervals (synchronous dy- 
namics). We specialize to the case where the positions Ui are integer- valued 
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and move by one lattice unit at a time: 



""'^^ ^ - \ n,{t) if /.(t)<0, 

where fi is the total force on site i. In the infinite-range limit it is given by 

fi{t) = u{t) - Ui{t) + Fext + g ViiUi), (9) 

the coupling parameter g controls the strength of the substrate disorder, the 
distribution p{r]) being normalized and of zero mean. 

Even for such a simple system it is not straightforward to determine the 
threshold very accurately, as for finite N the interface only travels a finite 
distance hm before it stops. This distance increases rapidly with in the 
mobile phase but it fluctuates strongly from sample to sample close to the 
threshold and the interface velocity is not precisely defined. The stopping 
distance plays the same role here as the polarization for charge-density waves 
p2[ . For a fixed driving force its average < hm > is expected to follow a 
scaling law of the form 

<hm>-^ H[Ny{g-g,)], (10) 

where the scaling function H{z) is finite and regular for z = 0, gdFext) is 
the critical coupling, and k and y are critical exponents. Expression ( [lO| ) 
contains three unknown parameters and it is more practical and accurate to 
study the Binder ratio 

Rh = <hl^> / <h^ >2, (11) 

which is expected to converge to a finite value at the critical point when 
N ^ oo. The results are displayed in figure (^ for a symmetric uniform 
disorder distribution, p{r]) = 1/2 for — 1 < r/ < 1. The critical ratio is 
R'^ ~ 1.28, a little below the value 4/3 which would correspond to a 
Gaussian distribution for h^, and this approach yields an estimate gc = 
2.380 ± 0.005 for the threshold at F^xt = 0. 

3.2 Exact solution 

An exact solution of the model may in fact be obtained by noting that deter- 
ministic evolution equations can be written down for Pk{x,t), the fraction 
of interface sites at height k which experience a local pinning force gx |21]. 
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Figure 2: Binder ratio Rh for the distribution of stopping distances as 
a function of the disorder strength g (1000 samples for sizes = 1000 to 
10000 and 100 samples of size 20000). 



The key is to realize that due to the simple dynamics this quantity has 
a quasi-factorized form : 

Pk{x,t) = l ^! ^<^^S' (12) 

^ ' ' \ l^k{t) p{x) if x>r]k{t), ^ ' 

where the discontinuity point is related to the mean interface position by 
= {k — u{t — l))/g- The weights Xk{t) and fJ-k{t) obey (relatively) 
simple recursion relations and the end result is a closed system of equations 
describing the interface dynamics in the thermodynamic limit. 

In the mobile phase {g < gc), this system admits only time-dependent 
solutions so the threshold may be determined by studying the existence 
of stationary solutions. For a uniform p{ri) the calculation can be done 
analytically and the threshold is given by the value of g for which a certain 
polynomial has a double root. For small driving forces one obtains pl| 



gc = 2.38006232... + 2.3901... F^^j. (13) 

Note that a threshold exists even for a negative Fg^t ■ Some interface ele- 
ments will move forward against the external bias, if the local pinning force 
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g r]i acting on them is strong enough, and due to the asymmetry of the dy- 
namics @ this is not compensated by backward motions of other elements. 
This "ratchet" effect is of collective origin - an individual particle in a neg- 
atively biased random potential would eventually stop -, and it is driven by 
quenched disorder rather than by random diffusion as in models recently 
proposed for molecular motors ||23[| . 

The threshold may also be obtained more directly through the following 
argument. Let us assume that in the pinning phase, among all the possi- 
ble equilibrium positions for an interface element, the one actually reached 
dynamically is the first lattice site z* for which /, is non-negative j^], i.e.: 

f{z*) = -z-z*+g^{z*) <0, (14) 
/(z) = z-z + gri{z) > 0, for ah z < z* , (15) 

where the position z of the center of mass can be restricted to < z < 1 
due to statistical translational invariance. The probability density of the 
stopping point is given by 

n^l = [\{ R{z^)][l-R{z*)l (16) 

Zi<Z* 

where R{z) is the probability that {z — z) < g r}{z). For a uniform symmetric 
distribution p{r]) one has 

{1 if z < z — 5, 

if z>z + 5, (17) 

i (1 - ^) otherwise. 

The self-consistency condition on the center of mass, z = z*V{z*), then 
yields a polynomial equation in z whose coefficients are functions of g and the 
condition for existence of acceptable solutions can be shown to be identical 
to the one obtained above for the threshold gc- 

The exact solution also shows that the interface has a bounded width 
and the average stopping distance < hm > is finite at g^ so 

C = K = 0, (18) 

in agreement with standard mean-field predictions. An analysis of the nu- 
merical results shown in figure (Q) gives a rough estimate y ~ 0.5 for the 
size exponent in (|lO|). 

Once the threshold is known exactly it is possible to study the inter- 
face dynamics in its vicinity in detail, analytically and numerically. One 
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finds that in contrast with the assumption usually made the velocity v{t) is 
strongly non-uniform: the interface moves through rapid bursts, separating 
periods of very slow motion. The minimum velocity Vmin vanishes linearly 
but the average velocity v, measured over a very long time interval, vanishes 
at the threshold as 

- 0.06811... (5c (19) 
which corresponds to a critical exponent 

6* = 1/2. (20) 

This result is surprising at first sight, as it differs from the standard 
mean- field prediction, 9mf = 1- It means that the existence of an underlying 
lattice is a relevant perturbation, in the RG sense. This can be understood 
physically, by noting that mean-field theory predicts the interface to be 
smooth (C = 0), i.e., its width remains comparable to the lattice spacing, so 
the continuum limit cannot be taken naively. In RG terms the discreteness 
prevents the effective random potential seen by the interface on large scales 
from developing a cusp at the origin . 

It remains to be seen if lattice effects also play a role for systems with 
long-range interactions such as the contact line. Numerically one finds that 
the crossover region between a linear and a square-root dependence of the 



form (19) is rather narrow, so it would be difficult to extract the correct 
critical behaviour from simulations such as those in section ( |3.1| ) - this should 
be kept in mind when analyzing numerical results in finite dimensions. 

I thank B. Derrida and A. Rosso for stimulating discussions, and M. 
Mezard for suggesting the alternative approach to derive the threshold value. 
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